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Abstract 

We analyse a binary mixture of colloidal parallel hard cylindrical particles with identical diam- 
eters but dissimilar lengths Li and L2, with s = L2/L1 = 3, confined by two parallel hard walls in 
a planar slit-pore geometry, using a fundamental-measure density functional theory. This model 
presents (|arXiv: 1002.0612V . accepted in Phys. Rev. E) nematic (N) and two types of smectic (S) 
phases, with first- and second-order N-S bulk transitions and S-S demixing, and surface behaviour 
at a single hard wall which includes complete wetting by the S phase mediated (or not) by an 
infinite number of surface-induced layering (SIL) transitions. In the present paper the effects of 
confinement on this model colloidal fluid mixture are studied. Confinement brings about pro- 
found changes in the phase diagram, resulting from competition between the three relevant length 
scales: pore width /i, smectic period d and length ratio s. Four main effects are identified: (i) 
Second-order bulk N-S transitions are suppressed, (ii) Demixing transitions are weakly affected, 
with small shifts in the /ii — /i2 (chemical potentials) plane, (iii) Confinement-induced layering 
(CIL) transitions occurring in the two confined one-component fluids in some cases merge with the 
demixing transition, (iv) Surface-induced layering (SIL) transitions occurring at a single surface as 
coexistence conditions are approached are also shifted in the confined fluid. Trends with pore size 
are analysed by means of complete //i — ^2 and p — x (pressure- mean pore composition) phase dia- 
grams for particular values of pore size. This work, which is the first one to address the behaviour 
of liquid-crystalline mixtures under confinement, could be relevant as a first step to understand 
self-assembling properties of mixtures of metallic nanoparticles under external fields in restricted 
geometry. 
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I. INTRODUCTION 

Understanding the self- assembly of nanowires and metallic or semiconducting nanorods 
will be crucial in next-generation nanoelectronic and display-technology applications [l|. 
These technologies exploit the ordering of nanorods to create well-ordered arrays for the 
selective fabrication of devices. Inorganic nanoparticles with controlled size, shape and 
composition, exhibiting interesting optical, magnetic and electric properties are now being 
synthetised 2|, |3| . But most applications are based on particle arrangements in large struc- 
tures. Confinement of the rods in nanoscale geometries may be one possible strategy to 
create such structures. External electric or magnetic field can be used to orient particles 
in the required directions [^. Predicting the shape and composition effects that lead to 
ordering in collections of nanometer-sized particles in a solvent is a theoretical challenge, 
and understanding entropic effects associated with the rodlike shapes is a prerequisite before 
considering other interparticle forces. Therefore, the interest in hard-particle models and 
in the statistical mechanics of positional and orientational ordering in fluids of identical or 
polydisperse hard particles has been revitalised in recent years. 

The conceptual aspects of freezing of monodisperse particles in nanopores has been stud- 
ied quite extensively J5l-ll6| . On the other hand, smectic and layered phases made of aniso- 



metric particles present partial (one-dimensional) ordering and give an opportunity to study 
partial, as opposed to complete, spatial order in conflned geometries [l7|. Frustration effects 
in mesophases, normally associated with antagonistic conditions imposed on the nematic 
director in the form of surface and bulk flelds, are crucial in technological applications of 



these materials [18| . But in smectic materials another type of frustration occurs, since the 
natural layer spacing d establishes a length scale which will compete with any other length 
scale applied in the form of geometric constraints. The simplest conflnement geometry is the 
planar slit pore, where the material is sandwiched between two parallel planar substrates. 
The conflnement of two- or three-dimensional one-component fluids in the smectic or lamel- 



lar regime leads to complex behaviour 19l-l21| induced by commensuration effects between 



the pore width h and the smectic spacing d. Theoretical work indicates that the possible 
phase behaviour includes conflnement-induced layering (CIL) transitions, as well as a 'wavy' 
nematic-smectic transition strongly coupled to the latter. The strong smectic regime has 
not been studied experimentally yet. 



In contrast to the case of pure materials, confined mixtures of liquid-crystal forming 



particles have been given no attention. In our previous work 2^, we have analysed the bulk 
and adsorption properties of a mixture of hard cylinders. A complex bulk phase diagram 
resulted, with nematic (N) and smectic (S) phases. The smectic phases came in two versions. 
In the first, non-microsegregated, smectic phase (Si), layers are identical and composed of a 
mixture, in various proportions, of the two species. In the second, or microsegregated smectic 
phase (S2), layers rich in one species alternate with layers rich in the other. Various phase 
transitions, including first- and second-order N-S transitions, as well as smectic demixing, 
were obtained; see Fig. [2I^a), where the bulk phase diagram is depicted in the plane pressure p 

lort particles). The adsorption properties 



22|. We found a complex behaviour, with 



vs. composition x (defined as mole fraction of the s : 
of the mixture on a hard wall were also studied in 
the wall inducing strong layering leading to stratification of the material into alternating 
layers rich in one of the components. The surface behaviour includes complete wetting 
by the S phase, mediated (or not) by an infinite number of surface- induced layering (SIL) 
transitions, as well as critical adsorption on approaching second-order N-S transitions. As 
interactions between particles and wall are of hard-core type, layering and wetting properties 
of the nematic fluid in contact with the wall have an entropic origin. Similar phenomena 
have been found in theoretical work describing the surface phase behaviour of colloid-polymer 



mixtures 23 | 



In the present paper we extend the analysis of this mixture by considering fluids confined 
by two identical hard walls. This is a novel system where a new length scale appears: apart 
from the pore width /i, there are two particle lengths, Li and L2. Commensuration between 
the three length scales brings about new phenomena. The most salient results of our work 
are: (i) Second-order bulk N-S transitions are suppressed, (ii) Demixing transitions are 
not much affected by confinement, with small shifts in the /ii — ^2 (chemical potentials) or 
p — X (pressure vs. mean pore composition) planes, (iii) CIL transitions occurring in the 
one-component fluids are modified by mixing, ending in critical points or at the demixing 
transition, (iv) SIL transitions on each wall are also modified by the confinement; the very 
first transitions (the lowest-order ones, involving the first few layers) may survive in wide 
pores, while the others are preempted by demixing; in some cases some of them may reappear 
in the demixing region creating islands of stability, while higher-order SIL transitions are 
suppressed. In general, both types of layering transitions tend to disappear as the pore 



width is decreased, (v) Some new phenomena, genuinely connected to confinement of the 
mixture, i.e. to the competition between the three length scales, Li, L2 and /i, may also 
arise. 

After summarising the density-functional model used in the analysis and discussing briefiy 
the bulk and surface phase behaviour in the next section, we present in Section UTTJ the results 
for the confined mixture, considering wide and narrow slit pores. Finally, some conclusions 
are presented in Section IIVI 



II. MODEL, THEORY AND BULK AND SURFACE PHASE DIAGRAMS 



A. Model 

We consider a binary mixture of parallel hard cylinders, with lengths Li, L2 and identical 
diameters D; see Fig. [U^a). Label 1 refers to the species of the short cylinders. We have 
selected Li as the unit of length. The length ratio investigated is s = L2/L1 = 3. The 
diameter D is adjusted so that the transverse particle area in units of Li is set to unity, 
i.e. TiD'^/4:L\ = 1. Since particles are parallel, bulk properties only depend on s, and 



the particular values of Li (or L2), and D are irrelevant. As discussed in J22|, the values 
of particle lengths chosen gives a binary mixture that would be equivalent to a mixture 
of freely-rotating but highly oriented cylinders (of the same length ratio) such that both 
components would have a stable smectic phase. The cylinders are in numbers A^i and N2, 
so that the composition of the mixture is defined as x = Ni/N, with N = Ni + N2 the 
total number of particles. The cylinder axes are chosen to lie along the z direction; see Fig. 
[U^b). This configuration models a real colloidal or molecular fiuid where, either by surface 
treatment or by means of a bulk field, particles are forced to point along some fixed direction 
which, in the case of an adsorption system, would be the surface normal. 

B. Theory 

The properties of the mixture are analysed using a fundamental-measure version of 



density- functional theory 24,1251]. This theory is in very good agreement with Monte Carlo 
simulations of monodisperse (one-component) cylinders, and there are no reasons to believe 
that the theory is not equally accurate for mixtures. We choose the grand potential Q as the 
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FIG. 1: (a) Schematic of particles studied, (b) Geometry of the adsorbed and confined fluid. 



relevant thermodynamic potential. If Pi{r) is the local number density distribution of the 
cylinder centres of mass of the i-th species, the grand-potential functional of the confined 
mixture is 

2 






(1) 



where /Xj is the chemical potential of the z-th species. The Helmholtz free-energy functional 
J-'[{pj}] is written, as usual [22|, as 



T[{p^}] = J^id[{pi}] + J^cA{P^}] + -^cxt[{Pi}]- 



(2) 



J-" consists of three parts: the ideal-gas part 
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the excess part. 



•FidKp.}] = kTj2 [ drPiir) {log [VMr)] " 1} 

-Fex[{p.J] = kT [ rfr<l>({p,(r)}), 
Jv 



(3) 



(4) 



and the external contribution from the walls: 

2 



-Fext[{p.}] = kTj2 f drp,{r)v}^{r). 



(5) 



Note that temperature T is not a relevant variable since all interactions are purely hard and 
therefore J-" is proportional to T. In the above k is Boltzmann constant, Vj is the thermal 
volume of species i, and $(|ft,(r)|) is the excess free-energy density. In the fundamental- 



measure approximation 



24 



25| . $(r) depends on the local averaged one-body nj(r) and two- 



body weighted densities A''j(r). However, by imposing the density profiles to be functions only 



on z (smectic symmetry), these weighted densities reduce to only two (one-body) densities, 
ni^z) and ?7(;z). In this case the expression for the free-energy density is given by 
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The averaged densities are 
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where cq = 7iD^/4 is the transverse area of the cyhnder. 77(2;) is the local packing fraction. 

The procedure to obtain the expression for $(r) consists of the following two steps: (i) a 
density functional is derived for a mixture of two-dimensional hard disks, with the important 
property of dimensional cross-over from two to one dimension (i.e. when disks are located on 
a straight line, the functional should reduce to the exact one for a mixture of hard segments), 
(ii) If we impose the parallel-alignment constraint on particles, the density functional for 
a mixture of parallel hard cylinders can be obtained from the density functional for the 
hard-disk mixture just derived by applying again the dimensional cross-over criterion. This 
is easy to visualize by taking into account that the projections of the cylinders on a plane 
perpendicular to their axes just give a mixture of disks with different diameters J24|. 

In the interfacial problems to be studied below, inhomogeneities resulting from the pres- 
ence of surfaces will be included via an external potential, which consists of hard walls: 



KiV) 



00, z < — or z > n , 

2 2 ' 



0, 



^ = 1,2, (8) 

2 - - 2 

where h is the width of the slit pore. The conditions on the confined fluid are controlled by 
the chemical potentials /ii and /i2 of the bulk phase with which it is in equilibrium. Equiva- 
lently the bulk phase can be characterised by the pressure p and the bulk composition x (in 
case of demixing in the bulk phase the transformation [fii, ^12) -^ {p,x) is not unique with 
respect to the composition variable). The bulk fluid mixture and the adsorption problem on 
a single hard wall (semiinfinite system) can be studied within the same scheme 22[. From 
the density profiles one can define the partial mean densities pi, the total mean density p 
and the global mole fraction x, as: 



Pi 



1 r'^ 

- / dzpi{z), 



P = pi +P2, 
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(9) 



In Sec. mil we will present phase diagrams of the confined mixture in the fii-fi2 plane but 
also in the p-x plane, since in the latter demixing regions can be appreciated. For fixed 
/ii, fi2, as /i — !> oo we have x -^ x. The equilibrium density profiles can also be obtained 
in other ensembles, and we have found it convenient in some cases to use the Gibbs free 
energy per particle G^ = G/N which depends, in particular, on the mean pore composition 
X. Gn allows the study of fractionation effects (gaps in x) at first-order transitions. This 
can be carried out by minimization of Gn with respect to p and Xi{z) = pi{z)/p for a fixed 
value of mean pore composition x. The Gibbs free energy per particle can be obtained 
from the Helmholtz free energy per unit volume Ty = J^/V by a Legendre transformation, 
Gn = (J-'y + r) /p, where T = — — — = — r2[{p,^'^''^}]/y, i.e. minus the grand potential 



dp 
per unit volume evaluated at the equilibrium density profiles [26[. Thus we obtain the usual 

result r = p, the bulk pressure, when external potentials are absent. The double-tangent 

construction on the function Gn{x) for a fixed value of F is equivalent to the equality of 

chemical potentials of different species in each of the coexisting phases. These chemical 

potentials in turn coincide with those appearing in the definition of f2[p]. Changing the 

value of F, and repeating the above procedure, we can calculate the phase diagram in the 

coordinates p-x. Note that the bulk pressure p can be calculated from the values of {pi}. 

In the following we summarise the bulk behaviour and the surface behaviour of the fiuid 

against a single wall. 

C. Bulk phase behaviour 

In 22| we analysed the bulk phases of the above model mixture. Here we briefiy discuss 
the main results. For the nematic phase pi{z) = pi =const., whereas in the smectic phase 
Pi{z) are periodic functions of 2, with period d. The stability of the nematic and smectic 



phases was obtained in a p-x (pressure-composition) phase diagram [22] (details of these 
calculations can be found therein). Fig. [2](a) presents the bulk phase diagram in this plane; 
in part (b) the same diagram is shown in the pi-p2 plane. The following features are worth 
mentioning: (i) there are two second-order N-S transitions (dashed lines), each starting at 
one of the pure-fiuid cases, x = or x = 1; (ii) one of these transitions, the one coming from 
the a; = axis, is connected to a first-order N-S transition by a tricritical point [triangle 
in Fig. [2](a)]; and (iii) three regions of S-S demixing, two of them ending in corresponding 
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critical points (circles), appear in the phase diagram. Note that, in the regions of smectic 
stability, two smectic structures occur: one where the density waves of the two species are 
in phase (one- layer smectic), denoted by Si, and another where the densities are out of 
phase (two-layer smectic), S2. The latter presents a microsegregation of the two species. By 
changing the conditions on the mixture one can pass from one of these smectic phases to the 
other in a continuous fashion. In the bulk phase diagrams of Figs. [2](a) and (b) , and in the 
rest of the article, the smectic phases are denoted by primed or unprimed labels according 
to whether the smectic is rich or poor in the short component, respectively. 

We now briefly comment on the expected changes in phase diagram topology if particles 
axes were allowed to rotate. The most important difference will be the presence of the 
isotropic phase, and the existence of an associated first-order isotropic-nematic phase tran- 
sition with a transition gap that depends on the mixture composition. Obviously our model, 
with the perfect-alignment constraint, cannot account for this transition. Also, our model 
predicts a second-order nematic-smectic transition which, as show below, is suppressed by 
confinement (an effect ultimately arising from the perfect-alignment approximation). How- 
ever, in situations where the bulk phases (nematic or smectic) are highly oriented, we expect 
a phase behavior very similar to that shown in the present work. 

D. Surface behaviour 

The behaviour of the single-wall system is very important to understand the properties 
of the confined mixture, as will be seen in the following section. The surface behaviour of 
the fluid when a nematic phase is adsorbed on a substrate was investigated in our previous 



work [22|. Particle axes were taken to point along the normal to the substrate, which was 
represented by an external potential Kxt(^) of the same type as in ([8]). We only studied 
the case where the bulk fluid in equilibrium with the surface structure was nematic. The 
nematic phase far from the wall {z — )■ 00) was characterised by values of the pressure p 
and the composition x or, equivalently, by the chemical potentials of the two species, /ii 
and p,'),. A conjugate-gradient scheme was used in the numerical minimisation (see details in 



Ref. [221]), from which equilibrium density distributions and values of surface tensions were 
obtained. 

At a given value of pressure p, wetting of the wall by the smectic phase was found for 
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FIG. 2: (a) Bulk phase diagram of the mixture in the "p-x plane. Text indicates type of wetting 
behaviour, (b) Bulk phase diagram of the mixture in the /ii-/i2 plane, (c) Phase diagram of the 
confined mixture in the p-x plane for pore width h = 48.2Li . (d) Zoom of the high p and x region in 
panel (c). (e) Phase diagram of the confined mixture in the p-x plane for pore width h = 6.8L1. (f) 
Phase diagram of the confined mixture in the fii-fj.2 plane for pore width h = 6.8L1. In all diagrams 
circles indicate critical points, triangles tricritical points, squares critical end points; dashed curves 
are continuous phase transitions. Dotted horizontal lines indicate triple or quadruple points. 
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all pressures as x — )■ x^s{p), where Xns(p) is the boundary of the N-S transition. However, 
depending on the bulk pressure, three regimes may occur: (i) Complete wetting by the 
Si phase, rich in long particles, via an infinite sequence of surface-induced layering (SIL) 
transitions when p is above the S1-S2-N triple point; at a first-order SIL transition, the 
composition of a localised interfacial region, with a thickness of one layer, changes abruptly 
from a low to a high value, as a result of which the short component is expelled from the 
wall and sharply localised, smectic-like layers, adsorb at the wall in a stepwise fashion, 
until a macroscopically thick smectic Si film develops at the wall (complete wetting), (ii) 
Complete wetting by the S2 phase without SIL transitions when p is below the triple point 
but above the tricritical point. And (iii) critical adsorption at the second-order bulk N-S2 
transition, i.e. for pressures below the tricritical point. The first of the infinite sequence 
of SIL transitions in the first regime has been plotted in Fig. Et^a); the rest are too close 
to the coexistence boundary to be seen. A blow up of the first four transitions is shown in 
Fig. E] in two different planes. Here, and in the rest of the paper, phases denoted by N„ 
refer to structures with n localised layers of the long particles near the wall, and a density 
that tends to a constant value far from the wall; in the confined case the density may not 
be completely uniform in the middle of the pore, but the structures are generally connected 
with the bulk nematic as /i — t- 00. 

III. CONFINEMENT 

One of the most striking phenomena presented by the mixture of cylinders is the occur- 
rence of SIL transitions at a single wall as the bulk N-S coexistence is approached. Our 
explanation for this phenomenon J22| is that there is an entropic competition (reduction of 
excluded volume) between the two species to cover the wall, which couples to translational 
entropy to produce a first-order phase transition. But there is another type of layering phe- 
nomena, i.e. the confinement-induced layering (CIL) transitions, which take place when a 
spatially ordered phase competes with a confining length scale. These phenomena are well 
documented and have been investigated theoretically in liquid crystals and related materials 



19 



20|. CIL phenomena should also happen in fiuid mixtures that can stabilise a smectic 
phase in bulk. Since the smectic phase of our mixture comes in two varieties (Si and S2) 
and the bulk phase diagram is already quite complex, we can anticipate a fascinatingly rich 
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(a) 




FIG. 3: SIL transitions in the planes (a) /U2-/U1, and (b) p-x in the single- wah system. In (a), 
AyU2 is the chemical potential of species 2 with respect to the value of the bulk nematic-smectic 
transition. In (b), Xcocx is the nematic value at the transition for each value of pressure. 

phase diagram when the mixture is confined and the possible occurrence of new phases. 

As explained above, the particle model used in our calculations exhibits a second-order 
bulk N-S transition. Due to suppression of fluctuations along the z direction when the 
fluid is confined, the N-S transition disappears in the confined system. Therefore, it is not 
possible to assign in a clear manner the roles of nematic and smectic to any of the confined 
phases, except in some particular cases (generally speaking, as pressure or chemical potential 
is increased, the oscillating distribution of particles tends to become sharper throughout the 
pore). The particular cases are: (i) Confined phases close to x > and a; < 1 existing in 
a closed region of stability, which are generally connected with the bulk smectic phase as 
/i — 7- 00. (ii) Surface phases in the confined case, which undergo SIL transitions, and which 
are connected with N„-type phases in the single-wall [h = 00) case (in bulk, these phases 
are in contact with a bulk nematic phase, but in the confined case the latter may have more 
or less spatial order depending on the values of pore width and pressure). Our convention 
is to denote all confined phases by S*-'"^ with m the number of layers of either species inside 
the pore, in case (i), and by N„, with n the number of absorbed layers of the long particles. 
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FIG. 4: (a) Gibbs free-energy density g* = G^ jkT as a function of global mole fraction x for 
the case /3pLf = 1.7. For the sake of visualisation, a straight line —10.96 + 6.OO2; has been 
subtracted from the free energies. Curves correspond to minima at fixed x obtained from the 
conjugate-gradient method (depending on the initial conditions different branches can be obtained). 
Continuous curves: branches that eventually lead to stable states. Dashed curves: metastable 
branches. Circles correspond to states obtained by simulated annealing. Dotted (straight) line 
indicates three-phase coexistence of a confined SIL transition, (b) Profiles of the metastable state 
indicated by the square in panel (a); the continuous curve is the profile of the long particles, 
whereas the dashed curve corresponds to the short particles. 

in case (ii). Intermediate cases are not well defined by these notations. 

In the confined case the method of solution of our theory is as follows. For fixed values 
of /ii, /i2, the fluid equilibrium structure within the pore was obtained using a conjugate- 
gradient (CG) method, where the z coordinate is discretised as z^ = kAz, with Az = O.OlLi, 
and p^ = Pi{zk) are used as discrete variables. However, the free-energy landscape of the 
confined mixture is very complex, with many local minima of similar free energy, correspond- 
ing to different arrangements of particles that can accommodate into the pore. The grand- 
potential minimum given by the CG method depends crucially on the initial configuration 
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FIG. 5: From (a) to (c), sequence of profiles in a simulated-annealing simulation for the case PpLf = 
1.7, X = 0.7 and h = 48.2Li. The continuous curves are the profiles of the long particles, whereas 
the dashed curves correspond to the short particles. A complete animation of the simulation can 
be found in 12711. 



used to start the iterative CG algorithm. Alternatively, we have used a simulated-annealing 
(SA) method to double-check that the minima found are actually the absolute minima. Even 
then, one always has to bear in mind that one cannot be completely certain as to the absolute 
nature of the minima found, especially in particular areas of the phase diagram where the 
particle arrangements may be more degenerate. This is the case in the region of the phase 
diagram where the values of p and x are high, i.e. in a compressed mixture where most of 
the particles are short. Fig. |4](a) shows the structure of minima in the Gibbs free-energy 
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landscape as a function of x for the case fipL\ = 1.7. Curves are minima obtained using the 
CG method, starting from different initial profiles. Depending on the initial condition, the 
minimisation procedure gets trapped into one of many possible metastable branches. Panel 
(b) of the figure shows a typical metastable state [actually the one located by the square in 
panel (a)]. For a given value of x, the different metastable branches correspond to the very 
many different ways of distributing the layers of the long species. Circles correspond to SA 
simulations starting from uniform profiles which hopefully, and in the cases explored it was 
always so, are able to identify the correct branch. In the case shown, there exists coexistence 
between three structures with different values of x, which follows from the double-tangent 
construction, for a SIL transition. Fig. \5\ shows the evolution of a SA simulation for the 
case (3pLl = 1.7, x = 0.7 and h = 48.2Li. The SA procedure is able to find the minimum 



free-energy configuration in a few hundred steps [27[. The resulting structure is then re- 
fined using the CG algorithm. Although the difference in free energy between the coexisting 
phases may be relatively small (as is the case shown in Fig. Hj), the values of the composition 
variable x are not so similar (see values at coexistence in the figure). This in turn means that 
the coexisting phases have different interfacial structures. Therefore, in case this scenario 
were to occur in an experimental system, large energy barriers separating these states may 
be expected, so that only fluctuations with relatively large amplitudes can take the system 
out from its original state. This scenario may be typical when large and moderate demixing 
gaps separate two stable confined phases. However, as the pore compositions of two coexist- 
ing phases become similar, and consequently also their interfacial structures become similar, 
these fluctuations, present in the real experimental system but not taken into account in the 
present model, can suppress these transitions. Inclusion of fluctuations in the orientations 
of the particle axes might stabilize interfacial structures such as those shown in Fig. Il](b), 
which is due to the lowering of the elastic energy when particles rotate in such a way as to 
modify the period of the confined smectic and better commensurate the number of smectic 
layers with the pore width. 

In order to understand the results for the confined mixture, let us first discuss the effects of 
confinement on the one-component fluid. The bulk fluid possesses a second-order transition 
from the nematic to the smectic, which disappears in the confined system due to the finite 
size available for spatial correlations. The only feature of the phase diagram that remains 
on confinement is the presence of first-order GIL transitions. These transitions occur at high 
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values of pressure or chemical potential, i.e. when the smectic order inside the pore is well 
developed. As mentioned already, their origin is the fact that, for a given value of pore width 
/i, the confined fluid can only accommodate a particular number n of smectic layers, which is 
given approximately by /i ~ nd, where d is the equilibrium layer spacing of the bulk smectic. 
In general the confined smectic will be stressed, either because it is compressed or swollen. 
As h is varied, the number of layers will be increased or decreased to optimise the available 
space from a thermodynamic point of view. Layering transitions are mostly associated with 
variations in /i, and depend more weakly on pressure or chemical potential, i.e. they are 
almost vertical in a p-h or ^-h phase diagram. CIL transitions terminate at critical points 



19|). There is a terminal, lower pore width at which layering transitions cease to 



(see Ref. 
exist. 



A. Wide slit pores 



Starting from the one-component fluid, as one adds a small amount of particles of the 
other component, the structure of CIL transitions is shifted in the phase diagram. If short 
particles are added to a smectic structure made of long particles, the smectic spacing will 
tend to increase and we expect the CIL transitions to occur at larger values of h. The 
distance between two consecutive transitions will be more or less constant. An example of 
this phenomenon is presented in Fig. El^a), which shows a CIL transition between smectic 
phases 83 and S-[ , containing 12 and 13 layers, respectively, in a mixture at low pressure 
{PpLf = 0.81). Note that there is a mole-fraction gap at the transition, i.e. the global mole 
fractions x of the two coexisting structures are different. There is also a terminal pore width 
beyond which the layering transition disappears. Above the transition gap in Fig. [6](a) the 
coexisting phase S^ smoothly changes to an out-of-phase structure 83 . On the other 
side of the phase diagram (x < 1), CIL transitions are not occurring at this value of the 
pressure since the mixture is in the N bulk regime, but they would certainly occur at higher 
pressure. 

At higher pressures the phenomenology becomes much more complex, as CIL transi- 
tions compete with demixing. This was to be expected, as the bulk phase diagram shows 
that demixing occurs for pressures /3pL'^ > 1.0. An example is shown in Fig. MJo). For 
the pressure shown, f3pLf = 1.51, strong N-S demixing occurs at bulk [Fig. [2]^a)], and 
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FIG. 6: Phase diagrams of the confined mixture in the h-x plane for two values of the pressure, (a) 
PpLf = 0.81. (b) PpLf = 1.51. Labels indicate nature of phases. Horizontal dotted lines denote 
triple points. 



correspondingly the confined mixture also exhibits demixing in a similar mole-fraction in- 
terval. But with important differences, which will be commented on later. First note that, 
in the range x < 1, smectic CIL transitions of the same type as in the region x > for 
the lower pressure do exist. They are separated approximately by a pore- width increment 
Ah ^ 1.2Li ^ d, the bulk smectic spacing of the short species. Fig. [7] shows the density 
profiles of two phases that coexist at the Sg S'j^ CIL transition. 

In the region x > 0, CIL transitions are rapidly preempted by a demixing instability. One 
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FIG. 7: Density profiles of two phases that coexist at a CIL transition in the confined fluid with 
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ppLf = 1.51 and h = 46.0Li [see Fig. EJb)]. (a) 82^ ^ phase, with x = 0.9385. (b) S^^ ^ phase, 
with X = 0.9391. The thin vertical line divides the slit pore in two identical parts. Continuous 
curve: species 2 (long particles). Dashed curve: species 1 (short particles). 

interesting feature of the phase diagram of Fig. [6](b) is that, in the proximity of the CIL 
transition, an asymmetric structure, denoted by S^l^, is stabilised in a small island in the 
h-x plane; presumably, these islands exist associated with other CIL transitions, creating 
phases Sn^m (in this notation, n and m are the number of layers of species 1 and 2 of the 
asymmetric structure, respectively). In the case shown, the island interacts with the CIL 
and demixing transitions, creating triple points where three confined fiuids coexist. For the 
S]^ -S^ CIL transition, the asymmetric structure, S^l^, contains 7 smectic layers of long 
particles at one wall, separated from 6 smectic layers of long particles at the other wall, by a 
single localised layer of the (minority) short species. Note that other very similar structures 
with almost identical values of free energy but with different positions of the single peak of 
the short particles may exist, but structures with more than one layer of short particles, for 
larger values of x, are unstable with respect to demixing. The triple points are indicated 
by horizontal dotted lines in Fig. Et^b), while the density profiles of the three phases that 
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FIG. 8: Density profiles of the three phases that coexist at a value of pressure /3pLf = 1.51 and 
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for a pore width h = 46.3Li [see Fig. [6)[b)]. (a) S\ phase, with x = 0.012. (b) S]^ phase, with 
X = 0.019. (c) S^l^, with X = 0.053. Continuous curve: species 2 (long particles). Dashed curve: 
species 1 (short particles). Note that, in panels (a) and (b), the density of short particles is too 
small to be visible in the graphs. 



j(14) 



j(a) 



coexist at the intermediate triple point, involving the S\ , S^^'*'' and S^'*;j'3 phases, can be 
seen in Fig. |S1 

Leaving aside these islands of stability within the demixing region, we can see that this 
region consists of two segregation 'bands': a first, wide band, in the mole- fraction interval 
0.04 < X < 0.65, and a second, narrow band, in the interval 0.68 < x < 0.76. The first 
band corresponds to the shifted bulk demixing N-S transition, while the second is associated 
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FIG. 9: Density profiles of the three phases that coexist at a value of pressure PpLf = 1.51 and 
for a pore width h = 45.0Li [see Fig. EJ^b)]. (a) Niphase, with x = 0.688. (b) No,i phase, with 
X = 0.724. (c) No phase, with x = 0.759. Continuous curve: species 2 (long particles). Dashed 
curve: species 1 (short particles). 

with the SIL transitions occurring in the single-wall system, but here in the confined case. 
Inspection of Fig. [9] indicates that this is indeed the case. In the figure, density profiles that 
coexist at a pressure PpLf = 1.51 in a pore of width h = 46.3Li are shown [filled circles in 
Fig. EU^b)]. Profiles in panels (a) and (c) are symmetric and correspond to the phases with 
low and high values of mole fraction, respectively. The structural differences between the 
two profiles are identical to those occurring at a single-wall interfacial structure undergoing 
a SIL transition, described in Sec. IIIDi The mole-fraction gap at this SIL transition arising 
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in the confined fluid tends to decrease as h is decreased, until it eventually disappears at 
a critical point, not shown in Fig. Mjo). The transition is connected to a corresponding 
SIL transition at a single surface as /i — ;■ oo. Inside the instability region of the confined 
SIL transition an almost vertical, continuous line has been drawn, which corresponds to a 
peculiar structure that coexists with the other two [panel (b) in Fig. [S]. This structure, 
which has interfacial structures identical to that of (a) at the left wall and to that of (b) 
at the right wall, is asymmetric but has the same interfacial free energy as the other two 
(an equivalent structure is obtained by interchanging the right and left interfacial profiles). 
The existence of structure (c) is the consequence of this phenomenon being associated to 
the SIL transitions, which in the confined case, and for wide pores, occur at both surfaces 
independently without any interaction. As the pore width is reduced, interaction between 
the two surfaces structures would make the line inside the instability region in Fig. M^h) 
to have a finite width. This structure seems too elusive to be detected experimentally, 
since domains of the two coexisting interfacial structures would probably mix in different 
proportions at both surfaces. 

The p-x phase diagram of the confined mixture for a particular value of pore width, 
h = 48.2Li (in the regime of wide pores) is plotted in Fig. Et^c). The diagram contains 
the same features as previously described, but viewed in a different way. Also, it allows 
to visualise critical points associated with confined SIL transitions and some additional 
features. The diagram is to be compared with the corresponding bulk diagram. The wide 
demixing region is quite similar. But some differences are apparent. First, as mentioned 
before, the second-order bulk N-S transition is vanished. Second, the structure of CIL 
transitions appears in the regions x > and x < 1, and is seen to either interact with 
demixing or terminate in critical points. Since h is fixed, the number of transitions in the 
phase diagram is very limited (CIL transitions occur mainly with respect to pore width and 
are weakly dependent on pressure). In Fig. [10] the density profiles of three representative 

(12") 

structures associated with the CIL transitions are shown. These phases are labelled as S2 , 
S2 and S-[ . In the first two, the density profiles of the two species are out of phase, 
whereas in the third the minority (short) species adopts an in-phase configuration with two 
small, separated peaks about the main smectic layer; in this case the short particles go into 
the main layer, forming a wide two-in-one layer where short particles fiow more or less freely. 
An important feature in the phase diagram of Fig. [21(c) is that the remnant of some of 
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FIG. 10: Density profiles of three representative states labelled in Fig. W[c) as Sg ,82 and S]^ 



^(14) 



:(12) 



(13) 



for a pore width h = 48.2Li. (a) S^ -' phase: /3pLf = 0.6, x = 0.14; (b) S^ ^ phase: l3pLf = 1.0, 
X = 0.10; and (c) S[ phase: /3pLf = 1.2, x = 0.05. Continuous curves: species 2. Dashed lines: 
species 1. 

the (infinitely many) SIL transitions occurring in the single- wall case are clearly apparent. 
For the pore width shown, h = 48.2Li, only four of them survive. The first one, involving 
the phases Nq-Ni, exhibits a long segregation region ending in a critical point; this structure 
is reminiscent of the corresponding Nq-Ni transition in the single-wall system, and in the 
variable x the transition shows a gap. The second transition, N1-N2, is much shorter but 
exhibits the same characteristics. The third transition occurs right in the middle of the 
demixing region and, instead of being a separated segregation region, it appears as a very 
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FIG. 11: Density profiles of three states that coexist at PpLf = 2.0 with pore width h = 6.8Li [see 
Fig. m^e)]. (a) sj^^ phase, with x = 0.944; (b) S^ phase, with x = 0.674, and (c) sf^ phase, with 
X = 0.010. Continuous curves: species 2. Dashed curves: species 1. 

small stability island [N3 in Fig. El^c)]. Finally, the N3-N4 transition is very short, and 
again presents an associated critical point. As mentioned in the discussion on the h-x phase 
diagram. Fig. |6]^b), asymmetric phases appear in the central region of the mole- fraction gap 
in the case of the first two SIL transitions; these structures are indicated by the labels No,i 
and Ni,2 in Fig. EJ^c). 

Finally, we focus on the region at the upper-right corner of the phase diagram of Fig. E^c). 
This region is quite complex. A zoom of this region is presented in Fig. Et^d). The small 
demixing region separating the bulk S'j^ and Sg phases, with the corresponding critical point, 
still survives for this pore width. Again, the bulk second-order N-S transition disappears, 
but the CIL transitions associated with the one-component short-particle fluid survive in 
mixtures with up to 12% of long particles. These transitions interact with the small demixing 
region. A final feature in Fig. [Sl^d) is the existence of a quadruple point, indicated by the 
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FIG. 12: Schematic of phase diagram evolution as the sht pore width is reduced, 
horizontal dotted line. 

B. Narrow slit pores 

In the case of narrow pores, available space is much more limited and much of the rich 
phenomenology described for wide pore vanishes. We have investigated the case h = 6.8L1. 
The p-x phase diagram is depicted in Fig. [2t^e). A large demixing region persists, but 
the diagram is otherwise relatively featureless: both SIL and CIL transitions seem to have 
disappeared, and the only feature that remains is a very narrow stability region, deep in 
the demixing region, where a largely asymmetric structure is stabilised. As a result, a new 
triple point arises. This asymmetric structure might be the remnant of a SIL transition in 
wider pores. Fig. [11] shows the density profiles of the three phases that coexist at the triple 
point. The phase shown in panel (b) may be interpreted as a demixed phase where mixing 
entropy is compensated by an optimised surface entropy. Fig. Mi) is the corresponding 
phase diagram in the /ii-/i2 plane. The demixing transition is slightly shifted with respect 
to the corresponding bulk transition [see Fig. [2](b)]. 

C. Kelvin equation 

The macroscopic approach to capillary transitions is given by the Kelvin equation, which 
relates the degree of under- (or over-) saturation of the capillary phase transition, with 
respect to a reservoir, in say chemical potential, A/x, to the surface tension 7^^ of the 
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interface between the two coexisting phases a and /3, the contact angle and the pore 
width h (in the case of a wetting layer, the pore width has to be diminished by twice the 
thickness of the wetting layer vS). The approach is only valid asymptotically, in the regime 
where h — )■ oo. For a mixture, the Kelvin equation has been discussed by Evans and Marini 
Bettolo Marconi 28|. We can apply it to our particular capillary demixing system. Let a. 
denote the phase favoured by the walls, in our case the 5*1 phase, and /3 the phase present 
in the central region of the pore, in our case N. We should remind ourselves that the real 
coexisting pore structure, N„, consists of two 'smectic' films adsorbed on the walls with n 
layers each and a central, more or less uniform, nematic film (the value of n depends on p 
and h. In a very wide pore, n will be large and it makes sense to split the N„ structure 
into two adsorbed films and a central one; in a narrower pore, the smectic layers in N„ may 
be just a few in number, and such a division, and as a consequence Kelvin equation, is not 
valid). Then the shift in chemical potential A/i2 of the long species (more abundant in the 



a phase) with respect to the bulk value is given by ! 



l-io^fiB ^cosG 

^/^2 = -j^^ , (10) 

hp\ 



where the B coefficient is 



B = ( ^-^ ] - ( L^] = ^^l_^^_ (11) 



X J ^ \ X J a 2^a2J/3 

and where the composition values are given at bulk. Since we have wetting at bulk coex- 
istence, cosG = 1 and h ^ h ~ 2w. To check the validity of the equation in our present 
setup, we have applied the equation for a pressure PpLf = 1.25 and for the wider pore width 
used in this work, hL^^ = 48.2 (in the case of the narrower pore, hL^^ = 6.8, the whole 
procedure is meaningless, as there is no way to identify smectic and nematic films in the 
fluid structure, since the pore is too narrow). We have w ~ 2L2 = 6 [for (SpLf = 1.25 we 
have coexistence between the Si and N2 phases, see Fig. EJ^c)], Xa = 0.0588, X/3 = 0.5823, 
13-fapLl = 0.055451, and p["^Ll = 0.0103, which gives B'^ = 0.0654 and /3A/i2 = 0.019. 
This is to be compared with /3A/i2 = 0.029 from the density-functional calculation. The dis- 
crepancy is a huge 35%, and the reason is twofold: (i) the pore width h = 48.2Li only spans 
approximately 15 lengths of the long particle; this is too small for the Kelvin equation to be 
accurate, (ii) In confined layered phases, such as the smectic phase, there occur important 
commensuration effects between the intrinsic periodicity of the phase and the pore width, 
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and an elastic contribution, containing; the layer compressibility modulus, to the Kelvin 
equation must be included |l7|, |l9| . For wide pores the simple Kelvin equation would give 
predictions for the capillary phase transition that are roughly averages over the real values 
(which are oscillatory with respect to the pore width, with a period of approximately one 
particle length). 

IV. CONCLUSIONS 

In this paper we have discussed the phenomenology of a confined mixture of hard cylindri- 
cal particles oriented in the direction perpendicular to the confining hard walls. Competition 
between the pore width h and the two lengths of the particles, Li and L2, create a complex 
behaviour. Three regimes can be identified, see Fig. [121 The first is the single-wall case, 
h = 00 [Fig. [TWa)]. where there is a wide demixing region and SIL transitions induced by 
the wall-fluid interaction; these off-coexistence transitions are associated with a regime of 
complete wetting by a smectic phase when nematic conditions prevail at bulk. In the second 
regime, that of wide but finite pores, Li,L2 <^ h < 00 [Fig. [TW b)]. a new feature arises: 
CIL transitions, induced by confinement, when the composition of the mixture is close to 
zero or unity (i.e. almost one-component cases), which may or may not interact with a wide 
demixing region. SIL transitions are affected by the confinement, and only a finite number 
of them survive. Finally, in the regime of narrow pores, i.e. wall separations only slightly 
larger than both particle lengths, /i > Li, L2 [Fig. [TW c)]. only demixing survives, and weak 
signatures of the layering transitions may be seen in the form of narrow stability islands in 
the demixing mole-fraction gap. As expected for mixtures of particles with such different 
volumes, demixing is always quite strong. The segregation region shows a surprisingly in- 
variant mole-fraction gap, and quite similar values for the pressure above which demixing 
occurs, with respect to pore width. This is reflected in Fig. [T31 where the demixing regions 
for the three cases investigated in this paper have been superimposed. 

This work has focused on confined phases with smectic symmetry. The possible stability 
of the columnar phase is a question that is left for future work. A complete calculation 
of the absolute stability of the columnar phase in the context of the present model is a 
difficult task which will have to be tackled using numerical techniques different from the 



ones developed here. In our previous work 



22| we used a bifurcation analysis to show that 
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FIG. 13: Demixing regions in the plane p-x for the three cases investigated in this work: h = oo 
(dashed curve), h = 6.8L1 (dotted curved) and h = 48.2Li (continuous curve). 

the nematic-columnar spinodal line is always below the neinatic-smectic line for all values 
of composition. This is not a proof that the smectic phases are stable with respect to the 
columnar phase, but at least indicates that part of the phases obtained here could be stable 
in some range of pressures. 
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